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Abstract. In this paper, we study an exact solution of the asymmetric simple 
exclusion process on a periodic lattice of finite sites with two typical updates, i.e., 
random and parallel. Then, we find that the explicit formulas for the partition 
function and the average velocity are expressed by the Gauss hypergeometric function. 
In order to obtain these results, we effectively exploit the recursion formula for the 
partition function for the zero-range process. The zero-range process corresponds 
to the asymmetric simple exclusion process if one chooses the relevant hop rates of 
particles, and the recursion gives the partition function, in principle, for any finite 
system size. Moreover, we reveal the asymptotic behaviour of the average velocity in 
the thermodynamic limit, expanding the formula as a series in system size. 



1. Introduction 

The asymmetric simple exclusion process (ASEP) is an exemplar of stochastic many- 
particle systems and provides fundamental models for various collective phenomena 
within the framework of nonequilibrium statistical mechanics [TJ |2j [3l 0]. With 
some particular conditions, the ASEP shows a wide variety of nonequilibrium critical 
phenomena such as boundary- induced phase transitions [5], spontaneous symmetry 
breaking [6], and phase separation [H El IS]- As described here, the ASEP has been 
extensively studied both with periodic boundaries and with open boundaries, and at 
the same time one may consider several update rules such as random and parallel [10J. 
In this paper, we focus on the ASEP on a lattice of finite sites with the periodic boundary 
condition and then consider a discrete-time evolution both with random dynamics and 
with parallel dynamics. The ASEP with the periodic boundary condition is suitable for 
applying a theoretical approach for the first time. 

In the ASEP, particles exclusively occupy sites on a lattice and hop from site to site 
in a definite direction with a constant probability, i.e., they move as if being driven by a 
biased field. The particles have a short-range interaction due to the exclusive occupation 
(i.e. the hard-core exclusion). However, as far as the parallel dynamics is concerned, 
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it may induce a long-range interaction as the density of particles grows. (Note that a 
parallel update rule treats all the sites equally, and the ASEP with parallel dynamics is 
hence regarded as a stochastic cellular automaton.) 

Generally in stochastic models, an asymmetric dynamics causes asymmetric 
transition rates between states of configuration and it is manifested as a current of 
probability [llj. According to that current, one observes a macroscopic phenomenon 
such as a density flow, and the system is hence thought to be far from equilibrium 
without any equilibrium state. However, nonequilibrium systems may have a steady 
state with a constant flow of particles, i.e., a nonequilibrium steady state [12]. The 
ASEP and the zero-range process, both of which are the subjects of the present paper, 
are typical examples for that kind of nonequilibrium system. We shall investigate the 
nonequilibrium steady state of these stochastic dynamical systems. 

Analytic methods were improved especially in the context of transport phenomena, 
e.g., traffic flow, and then exact solutions were presented by powerful theoretical methods 
such as the matrix-product ansatz [13], and the cluster-approximation method [14J. 
However, the previous results are for systems in the thermodynamic limit, i.e., of infinite 
system size. As far as nonequilibrium many-particle systems are concerned, it is essential 
that one should investigate these systems of finite size and then consider them in the 
thermodynamic limit. In this paper, we present an explicit formula of the partition 
function for the ASEP on a ring of finite sites, which allows for an exact calculation of 
physical values. 

We find that the partition function is expressed by the Gauss hypergeometric 
function. In order to obtain this result, we exploit the recursion formalism for the 
partition function for the zero-range process (ZRP) [15], which is given in [16]. As 
described later, the ASEP corresponds to the ZRP if one chooses the relevant hop rates 
of particles. Since the partition function for the ZRP is, in principle, calculated for any 
finite system size via the recursion, one is able to obtain that for the ASEP in the same 
way. 

As an application of the present result, we calculate an explicit formula for the 
average velocity of particles in the ASEP. Due to the formulas for the hypergeometric 
function, we can provide an explicit formula for the average velocity by the Gauss 
hypergeometric function as well as for the partition function. In addition to that, we 
provide the series expansion of the average velocity in system size by solving a sequence 
of equations of increasing order. The series expansion reveals the asymptotic behaviour 
of a physical value in the thermodynamic limit. In particular, we recover the well-known 
result given in [T4] . 

In most recent works [TTJ [18] , the same formulas as we obtain in the present work 
are given by using the Bethe ansatz solution of the master equation for the ZRP. They 
find that the hop rates required by the applicability of the Bethe ansatz allow two 
parameters and in a special case the ZRP reduces to the ASEP. It is interesting that, 
although the ZRP is concerned in both approaches, they seem entirely different until 
just before reaching the same formula. We remark that, in contrast with the Bethe 
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ansatz, our approach can be generalized to provide an exact solution of the ZRP if one 
takes other hop rates which do not comply with the conditions required by the Bethe 
ansatz [T9] . 

This paper is organized as follows. In section 2 we define the ZRP and see that in a 
special case the ZRP corresponds to the ASEP. Then, we derive the partition function 
for the ASEP after preparing the general formalism for that for the ZRP. In section 
3 we provide the series expansion of the average velocity in system size. Section 4 is 
devoted to conclusions and final remarks. Long derivations of formulas are included in 
the appendix. 

2. Exact solution of the asymmetric simple exclusion process 

2.1. The zero-range process 

The zero-range process is an exactly solvable stochastic process, and it has been 
widely used as a model for systems of many particles interacting through a short-range 
interaction [15J. Particles in the ZRP are indistinguishable, occupying sites on a lattice 
(any lattice in any dimension), and each lattice site may contain an integer number of 
particles. These particles hop to the next sites with a rate which depends on the number 
of particles at the departure site. One is to apply the ZRP to phenomenological studies 
choosing a relevant hop rate of particles which determines the microscopic behaviour. 

One of the most distinct properties of the ZRP is that its steady state is given 
by a factorized form [15]. Each of these factors of the steady state are for each site of 
the lattice, and are determined by the hop rates. In this sense, these factors are called 
the single-site weights. We remark that the steady state behaviour is determined by 
the single-site weights and consequently one may infer the hop rates after choosing the 
single-site weights as desired. 

The ZRP in one dimension can be mapped onto the ASEP if one rearranges particles 
and sites in the ZRP on a lattice in the ASEP in order that site m, containing n m 
particles in the ZRP, should be thought of as, in the ASEP, the mth particle occupying 
a site and unoccupied n m sites, i.e., the distance to the (m — l)th particle in front of 
the mth particle. In particular, the ZRP corresponds to the ASEP if one chooses the 
relevant hop rates (defined by ( fl4T) ). In recent works [201 EI], we proposed a traffic-flow 
model connecting the ASEP and the ZRP with an additional parameter. 

In what follows, we consider the ZRP in one dimension with N particles on a 
periodic lattice containing M sites labelled I — 1,2, ... , M. It corresponds to the ASEP 
with M particles on the periodic lattice of L (= M + N) sites, and then one should 
understand that the system is of size L and the density is p = M/L. 

2.2. The partition function for the ZRP 

In this subsection, following |15| . we define the partition function for the ZRP and then 
find that the generating function for the partition function is obtained from that for the 
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single-site weights. 

The nonequilibrium steady state probability P({n m }) of finding the system in a 
configuration {n m } = {ni,ri2, ■ ■ ■ ,%} is given as a product of the single-site weights 
denoted by /(n): 

^ M 

P({nm}) = -z II f( n m) (n 1 +n 2 + -- + n M = N) 1 (1) 

where a normalization Zm,n is the so-called partition function. The probability that a 
given site (e.g. site 1) contains n particles is given by 

p(n)= £ P({n,n 2 ,...,n M }) = f(n) ZM - 1 ' N - n . (2) 

ri2+n s + -+n M =N-n ^M,N 

The sum of p(n) over n is unity by definition. Thus, we obtain the recursion formula 
for the partition functions, 

N 

Zm,N = £ f( n )ZM~l,N-n, (3) 
n=0 

Z 1>k =f{k) (k>l). (4) 

Note that this is a recursion for a double series with respect to M and N, and the 
partition functions are calculated recursively from the initial values Z\ ^. 

Considering the generating functions /(£) := Y^=of{ n )C n ^ an d Z M (() : = 
J2^Lo ZM,n( n , w e have a recursion with respect only to M: Z M (() = f(()Z M _i((). 
Accordingly, we find the fundamental relation 

Zm(C) = (/(Of ■ (5) 
Thus, the partition function Zm,n for auy N is obtained from the single-site weights 
/(")■ 

The single-site weights, being obtained from the hop rates of particles, entirely 
characterize the zero-range process. In [T5l [T6] , two formulas of the single-site weights 
respectively corresponding to random and parallel update rules are given as follows. For 
a random update rule, the single-site weights f{n) are expressed as 

n i 

/(0) = 1, f(n) = l[—- (n>l), (6) 

where u(n) are the hop rates of particles when the departure site contains n particles. 
Note that u(0) = by definition. For a parallel update rule, the single-site weights are 
expressed as 

[ l-u(l) (n = 0) 

/(„) = 1^11 g 1^ > (7) 
[ 1 - u{n) u{j) 

It is remarkable that according to the update rule the single-site weights have 
different recursions. From one finds that the single-site weights for a random update 
rule have a recursion 

u(n)f(n) = f(n-l). (8) 
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In a similar fashion, from §J§ we have the recursion that the single-site weights for a 
parallel update rule satisfy, 

u(n + l)/(n + 1) = f{n) - u(n)f(n). (9) 

We remark that as for the single-site weights the recursions ([8]) and (Q may identify the 
different update rules instead of the explicit formulas ([6]) and ([7]), which will be seen in 
the course of calculations on the average velocity of particles in the ZRP. 



2.3. The average velocity of particles in the ZRP 

In this subsection, we shall provide the explicit formulas of the average velocity (or the 
mean hop rate averaged in the nonequilibrium steady state) of particles in the ZRP 
with random and parallel update rules. In general, the flux of a transport system briefly 
presents its macroscopic property, and especially in nonequilibrium systems it is one of 
the few quantitative criteria. In the present study, the flux of particles is obtained by 
multiplying the average velocity by the density of particles, since the number of particles 
is conserved due to the periodic boundary condition. 

The average velocity, denoted by vm,n, is defined by 

N N Z 

v M ,n = u ( n )p( n ) = H u ( n )f( n ) M ~ 1 ' N ~ n . (10) 

n=0 n=0 ^M,N 

As mentioned in subsection \2.2\ the recursion for the single-site weights plays the central 
role in the subsequent calculations. 

We first consider the ZRP with the random rule. From ([3]), ([8]) and ({TO]) , we find 
that the average velocity is expressed by the partition function, i.e., 

vm,n = • (11) 

Next, we consider the parallel update rule. From ([3]), ([9]) and ffTU]) . we find the 
average velocity satisfies the recursion with respect to N, 

vm,n+iZm,n+i = Z M ^ N — vm,nZm,n- (12) 

(It is suggestive that ( Ti~2l is corresponding to (JSJ) -) Solving (fl2i) for vm,n with respect 
to N, the average velocity is expressed by the partition function, i.e., 

Vm,n -, : s Nv • K*- 6 ) 

2.4- Exact solution of the ASEP 

Now, we turn to the calculations for the partition function and the average velocity in 
the ASEP. As mentioned in subsection 12.11 one chooses the hop rates in the ZRP as 

u(0)=0, u{n)=p (0 < p < 1, n > 1), (14) 

and the ZRP is thereby conformed to the ASEP with hop rate p. 



Exact solution and asymptotic behaviour of the asymmetric simple exclusion process on a ringQ 



J M,N 



17) 



2.4-1- Random update rule First, we consider the ASEP with the random update rule. 
From ([6]), the single-site weights are f(n) = p~ n (n > 0) and then the generating 
function for them becomes 

*>=sG)" = i^ (15) 

Hence, from <§5§ the generating function for the partition functions is figured out as 

Zu(0 = fi - C -) M = E ^r?C n (16) 

Thus, we obtain the partition function for the ASEP with the random update rule, 

(M) N 1 (M + N-l\ 
p N N\ ~ p^y M -1 J ' 

where (a) n = a(a + 1) ■ • • (a + n — 1) is the Pochhammer symbol. Therefore, from ffTTl) 
and f[T7|) the average velocity for the ASEP with the random update rule is given by, 

Np 

VM > N = M + N-l (18) 

2-4-2. Parallel update rule Next, we consider the ASEP with the parallel update rule. 
From ([5]) and ([7]), the generating function for the partition function becomes 

Expanding (fl9l) in a power series of (, we obtain the partition function for the ASEP 
with the parallel update rule, 

Zm» =(!-?-)", (20) 

^ j-fH^ /MH^.n (JV21) _ (21) 



where 



'a,/3 \ ^(a)„(/9)n^ 



2; 



E (22) 



7/^0 (7)n « 



is the Gauss hypergeometric function. (See Appendix A for details of the above 
calculations.) This presentation allows one to take advantage of the formulas for the 
hypergeometric functions to make advanced calculations. 

In order to calculate the average velocity according to (|T3l) . we first evaluate the 
sum in the numerator. Using the Gauss recursion formula with respect to parameters 

-M a '^^V4 (23) 



one finds 

AT-l 



n=0 \ 1 1 PJ 
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Finally, the average velocity for the ASEP with the parallel update rule is obtained as 



F 



(M,N l 



p — 1 y 1 ' 1 — p J 

" M - N = ~irjM+i,N+i i v (25) 

[ 2 '—pj 
3. Asymptotic behaviour of the exact solutions 

In this section, we investigate the asymptotic behaviour of the exact solutions of the 
ASEP, i.e., we find the power series expansion of the partition function Zm,n and the 
average velocity vm,n, for a given density p = M/L, with respect to the system size 
L (= M+N). As mentioned in the introduction, the average velocity is a physical value 
and is hence expanded in a power series with respect to the system size. In contrast, 
the partition function shall be expanded as an asymptotic series. 

3.1. Random update rule 

First, we consider the random update rule. After being expressed with the system size 
L, the average velocity ffTTj) is easily expanded in a power series as follows. 

v m ,n = P{1 ~_ P 1 L = P(l -/>)(! + L- 1 + L~ 2 + ■ ■ ■)• (26) 
From ( |26i) . we have an expected consequence 

lim v m ,n = pO- - p). (27) 

3.2. Parallel update 

Next, we consider the parallel update. For the sake of convenience, we shall present the 
power expansion of the average velocity before turning to the partition function. To 
begin with, we change the independent variable from p to z = 1/(1 — p). Then, using 
the formula for the derivative of the hypergeometric function: 
d n 

7 



dz T > 



1<-\1 - z ) a +P-yF^ a [P; z 



(28) 



= ( 1 -n) n z^- n -\l-z) a+ ^- n F( a n,P n ;z), (29) 

\ l-n J 

we can transform the rhs of ff25l) into a logarithmic derivative: 

z-l d , , „ /M + l, N + l 



VM,N 



M dz 



log (^(l — z) M+N F ( ' ; zn. (30) 



The Gauss hypergeometric differential equation for the hypergeometric function in the 
argument of the logarithm in fl30|) is given by 

z(l — z)—^- + [2 — (M + N + 3)z)^- - (M + 1)(N + l)w = 0. (31) 

dz z dz 
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Accordingly, the argument of the logarithm also satisfies a differential equation 

d 2 w + 1- M-Ndw + MN w _ q ^ 
dz 2 z — 1 dz z(z — 1) 

and we thus find that the average velocity vm,n is the solution of a Riccati equation 
P(p ~ 1) d 2 fi,H\ Np 

~^ Tp VM, N = v M>N - (1 + - j , M ,iv + W (33) 

Expanding the average velocity as vm,n = vo(p) +vi(p)L~ 1 + V2(p)L~ 2 + ■ • -, we separate 
the Riccati equation according to power of L: 

w o - ~ v o + = 0, (34) 

P P 

PiP ~ = Yl P v k y i ~ v i U - x )- ( 35 ) 

dP fc+i=j 

fc,i>0 

This sequence of equations, solved in turn starting from vq, give the power series 
expansion of the average velocity Vm n with respect to the system size L, 



Vm,N 



1 - V 1 - 4 PP( X ~ P) (l-p)p(l-p) ! 
2p 1 - 4pp(l - p) 

1 _ p ) p (i _ p ) [i _ 2 p + p(3p + l)p(l - p) 



+ 



l-4pp(l-p) 



+ •••• (36) 

In particular, we recover the well-known result for the average velocity in the 
thermodynamic limit [H], i.e., 



, 1- y/l-4Ml-p) _ 
hm w MjiV = 1 . (37) 

L^oo 2p 

4. Conclusion and remark 

In the present paper, we provide an exact solution of the ASEP on a periodic lattice 
with respect to two typical dynamics, i.e., random update and parallel update. To begin 
with, we focus on the following two facts; First, the ZRP corresponds to the ASEP if 
one takes the hop rates of particles to be a constant in the ZRP. Second, the partition 
function of the ZRP is exactly solvable in the sense that the partition function for any 
system size is obtained by a recursive calculation. Then, by solving the recursion with 
the relevant hop rates, we are to obtain the partition function for the ASEP. 

For the random dynamics, it is given by the product of a power of the hop rate and 
the binomial coefficient in which the number of sites and that of particles appear. For 
the parallel dynamics, it is given by the Gauss hypergeometric function in which the 
number of sites and that of particles equally appear as the parameters. It is remarkable 
that the hop rate becomes the independent variable. Thus, the mathematical structure 
of the ASEP becomes clear from the viewpoint of special functions. 
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By using the above results, we calculate the average velocities and obtain their 
series expansions in system size. The series expansions confirm the previous results 
given in the thermodynamic limit. Note that, if necessary, one can obtain the higher- 
order correction terms for the average veloctiy in the parallel dynamics by solving the 
sequence of equations in order. 

From the present results, it is expected that one can estimate the partition function 
for a ZRP with other hop rates than those we take in this paper. This will be reported 
in the near future. 
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Appendix A. Calculation on the generating function for the ASEP 

From the hop rates (JHj), the single-site weights are 



V P J 

For simplicity, we let 7 = 1 — p and (3 = (1 — p)/p hereafter. Then, the generating 
function for the single-site weights is 






/(o = 7+E(/?cr = 7+ T 

n=l 1 




(A.2) 



Accordingly, the generating function for the partition function becomes 




(A.3) 



Using the Euler transformation for a series: 




b 



11 



where 




(A.4) 



00 




(A.5) 



n=0 



n=l 



where 




(A.6) 
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From (IA.6p . we have 




r=0 



(2)r(l) 

After all, we obtain flT2J) from flA~5D and flA~9|) . 



(A.7) 

— L J \ V I 

me 1 "-:!' 1 '-^^-' (a.8) 



^MFl 1 - M - 1 - lf ;r l ). (A.9) 



Moreover, using the Rummer's transformation formula 

f( z\= (1 - zy-^F^ - a ^ -P-z\ (A.IO) 
one obtains (TT9D from (EH) . 
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